FLUX-LATTICE MELTING IN TWO-DIMENSIONAL DISORDERED 

SUPERCONDUCTORS 
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The flux line lattice melting transition in two-dimensional pure and disordered superconductors 
is studied by a Monte Carlo simulation using the lowest Landau level approximation and quasi- 
periodic boundary condition on a plane. The position of the melting line was determined from the 
diffraction pattern of the superconducting order parameter. In the clean case we confirmed the 
results from earlier studies which show the existence of a quasi-long range ordered vortex lattice 
at low temperatures. Adding frozen disorder to the system the melting transition line is shifted 
to slightly lower fields. The correlations of the order parameter for translational long range order 
of the vortex positions seem to decay slightly faster than a power law (in agreement with the 
theory of Carpentier and Le Doussal) although a simple power law decay cannot be excluded. The 
corresponding positional glass correlation function decays as a power law establishing the existence 
of a quasi-long range ordered positional glass formed by the vortices. The correlation function 
characterizing a phase coherent vortex glass decays however exponentially ruling out the possible 
existence of a phase coherent vortex glass phase. 

PACS numbers: 



-a 
o 

CN 
> 

CN 
CN 
CN 

CN 
O 
-i— > 



i 

C 

o 
o 



X 



I. INTRODUCTION 



The nature of the mixed phase in type-II two- 
dimensional (2D) superconductors, despite considerable 
experimental, theoretical and numerical effort, remains 
unclear. In the mean field approximation the vor- 
tices of a pure system are known to form an Abrikosov 
lattice^. In D=3 dimensions thermal fluctuations re- 
duce the upper critical field line as first shown by 
Eilenberger— leaving a vortex liquid phase between the 
transition line and the mean field H C 2(T). In D—2 di- 
mensions the phase diagram between the critical field 
lines H c i(T) and H C 2(T) is still under debate. Apply- 
ing the Kosterlitz-Thoulcss-Bcrczinskii-Halpcrin-Nelson- 
Young (KTBHNY) theory 3 i 4 i 5 i 6 i 7 originally developed for 
melting in 2D crystals to superconductors^ one ob- 
tains a phase diagram with a solid, a hexatic and a liq- 
uid phase. The hexatic phase was observed in recent 
simulations^ for the pinned XY model and low exter- 
nal magnetic fields. However, the possibilty that the 
transition is not of KTBHNY type was not excluded 9 -. 
Indeed, for higher densities of dislocations (i.e. their 
core energy is low) a first order transition was foundii. 
Moreover, vortices are not hard core particles. This ap- 
plies in particular if they are described in the lowest 
Landau level(LLL) approximation where the density of 
vortices is high. With the help of a high temperature 
series expansion starting with the LLL approximation 
Hikami et a! 12 found indeed a first order transition from 
the quasi-long range ordered vortex lattice to the liq- 
uid phase. In this case the intermediate hexatic phase 
disappears. The same result was obtained by Tesanovic 
and Xingi^ by Monte Carlo simulations using the LLL 



approximation 1 ^^. However, theoretical works of Moore 
and coworkersi&i! raise doubts about the existence of 
the quasi-long range ordered vortex lattice and suggest 
that the flux liquid phase exists at any nonzero temper- 
ature. The simulations based on the LLL approximation 
gave conflicting results regarding the flux lattice melt- 
ing in a clean 2D system. Using quasi-periodic boundary 
condition on a plane, Hu and MacDonalcU^, and Kato 
and Nagaosai 9 - have found a first order transition from a 
quasi-long range ordered vortex lattice phase to a vortex 
liquid. However, the quasi-solid phase was not detected 
by simulations on a sphere2fl*2ii22. 

In the presence of disorder the situation becomes even 
more complicated. It has long been predicted that in 
the presence of frozen-in disorder any long range crys- 
talline order of a vortex lattice phase is destroyed in less 
than four dimensions 2 ^. Instead, a liquid-like phase was 
expected to occur leaving only short range crystalline or- 
der at length scales smaller than the Larkin length 24 . 
Recently, it was suggested that a Bragg glass phase with 
quasi-long range order can exist in three dimensional im- 
pure superconductorsSSiSLSL 2 ^ provided the disorder is 
weak enough such that the dislocations cannot prolif- 
erate. Its existence in three dimensions was supported 
by further analytical22^2i2i and numerical^ 2 ^ studies. 
Experimental evidence for the existence of the Bragg 
glass phase was provided recently by neutron diffraction 
data2i. 

The situation in the 2D case is less clear. For clean sys- 
tems the KTBHNY theory of melting predicts a disloca- 
tion driven melting transition from a positionally quasi- 
long range ordered to a hexatic phase. The transition 



temperature is T m = 



47T A+2/j ' 
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and /j, = C66 denote the effective Lame coefficients, 

denotes the flux quantum. 
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For T < T m the structure factor 



5(k = G + q) 



|2-r, G (T) 



(1) 



shows quasi-Bragg peaks with a temperature dependent 
exponent 



Vg(T) 



|G| 2 TA + 3^ 
4tt^i A + 2[i ' 



(2) 



In the 2D Abrikosov triangular lattice where the screen- 
ing length diverges with vanishing thickness of the film, A 
becomes scale dependent X(L) ~ L 2 , implying r]G(T m ) = 

l G | 2 °A _ 1 QQ|2 _ 16tt 2 



-). In finite systems of size L s 



16tt 2 3 \l I 3a^ 

VG(Tm) may be larger because the effective value of X(L) 
is still finite. Above the melting transition the transla- 
tional long range order of the lattice decays on scales 
L > £ ss exp [ J_5? ] with i/ ss 0.3696 7 . In finite sys- 
tems one will see a slightly larger melting temperature 
defined by L s £,{T). The hexatic phase shows still 
quasi-long range order of the bond orientations which 
will disappear at a second transition at higher temper- 
atures. As mentioned already, the possibility that the 
transition from the crystalline to the liquid phase is first 
order cannot be excluded. 

Adding disorder to a purely elastic 2D system, i.e. if 
one suppresses dislocations by hand, there is at low tem- 
peratures a glassy phase with a structure factor 



S(k = G + q) 



|2-r, G (T)ln(|q|fl a ) 



(3) 
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where fjci(T) ~ _/ g V with T g being the glass tem- 
perature and R a is the length on which the vortex dis- 
placements become of the order of the lattice spacing 
a a- One should also note that |q| has to be smaller than 
L" 1 where L co is a cross-over length scale to the asymp- 
totic behavior. The phase characterized by the structure 
© was found by Carpentier and Le Doussal2£ following 
earlier work by Cardy and Ostlunoj 3 ^. We will call this 
phase therefore the Cardy-Ostlund-Carpentier-Le Dous- 
sal (COCD) phase. Contrary to the Bragg glass phase 
in D=3 dimensions this phase does not show infinitely 
sharp Bragg peaks for |q| — > but the Bragg peaks satu- 
rate at qR a < e~ 2 / flG as can be seen from Q. Moreover, 
this phase is not stable with respect to dislocations which 
appear on a length scale 3 ^* 3 ^ 3 ^ 



d[2 



■ ln(R/R a )] 1/2 
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due to the effect of the disorder. Here c\ is a numerical 
constant. Since in the thermodynamic limit the glass 
temperature T g — p^jgj~y is much higher than the 

melting temperature T m (y 2 -^ ^ x — 6) one will never 
observe the glass transition. Thus, for system sizes L s < 
Ldis we have to expect to see essentially a temperature 



driven melting transition from a low-T COCD phase to 
a vortex liquid. 

It should be noted that the LLL simulations 3 ^ for dis- 
ordered systems fail to see any evidence for the existence 
of the COCD phase in even weakly disordered 2D su- 
perconductors of spherical geometry. The same problem 
but for a plane with the quasi-periodic boundary condi- 
tion has not been studied yet. 

In this paper we employ the LLL approximation and 
quasi-periodic boundary condition on a plane to study 
disordered 2D type-II superconductors by Monte Carlo 
(MC) simulations using the standard Ginzburg-Landau 
model where the impurities are introduced via fluctua- 
tions of the local critical temperature. The disorder is 
assumed to be weak such that the distance L^is between 
dislocations is much larger than system sizes studied and 
hence the dislocations can be excluded. 

In order to study the glass phase we monitor the 
Fourier transform of the density-density correlation 
functioni&ia 



Xdd (k) = J J < |*(r)|2|*(r')| 2 >e^ r - r '\ 



(4) 



where ^(r) = l^le^ denotes the superconducting order 
parameter and the overbar the disorder averaging, re- 
spectively. The value of this function at k = G with 
G being a reciprocal lattice vector defines the intensity 
of Bragg peaks. A closely related correlation function is 
that of the order parameter for translational long range 
order defined as follows 37 



S(G,r) = <exp(iG.[u(r)-u(0)])>. 



(5) 



Here u(r) denotes the displacement field of the vortex 
positions. The latter follow from the condition ^(r) = 
0. It is this correlation function which exhibits quasi- 
long range order in the Bragg glass phase. The structure 
factor is the Fourier transform of S(G,r) 



d 2 re l 



r 5(G, 



In analogy with spin glass theory^ one may further 
consider the positional glass correlation function 37 



S PG (G,r) = I < exp(zG.[u(r) - u(0)]) > 



(6) 



which may give signature of the existence of some resid- 
ual 'glassy' order of the Abrikosov lattice. If 5pc(G,r) 
decays not faster than a power law for |r| — > oo then a 
system is said to be in the positional vortex glass phase. 

In analogy to the positional glass correlation function 
one can define the gauge-invariant phase-coherent vortex 
glass correlation functionii Cvg(?)'- 



Cvc{r) = | < *(r)**(0)e 4 ( 27r /*°)i'r £;r - A > | 5 



(7) 



Note that Cvc( r ) itself depends on the path T between r 
and along which the vector potential is integrated^. To 
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make this correlation function path independent and si- 
multaneously gauge invariant, Moored proposed to keep 
only the longitudinal component of the vector potential. 
This corresponds in the symmetric gauge to the complete 
neglect of the phase factor exp [j(27r/$ ) / r dr. A] in Eq. 
0: 



C VG (v) = | < ¥(r)¥*(0) 



> 



(8) 



The correlations measured in © cannot be probed di- 
rectly but have been proposed to characterize the glassy 
state of the vortex array. They are expected to be re- 
flected however in the current voltage relation although 
a straightforward derivation of this connection is still 
missing 41 . The relation between the correlations mea- 
sured by an d the dissipation due to flux creep is 
closer—. Since we do not discuss dynamical properties in 
this article we consider © and JSJ (and JBJ as two com- 
peting theoretical tools to characterize the glassy state of 
the vortex array. An exponential decay of © - as we find 
- would suggest , according to 41 -, the disappearence of the 
glassy state. This is in contrast to what one would con- 
clude from the almost power-law like decay of (J5J which 
suggest some remaining glassy COCD-behaviour in the 
low temperature state. Contrary to previous works we 
will focus on the behavior of correlation functions defined 
by Eqs. ©, © and ©. 

In agreement with reii^i^ we find that the clean 2D 
system displays the quasi-long range ordered vortex lat- 
tice phase at low temperatures. This result is in variance 
with the result of Moore et a&SLS^ who found the vor- 
tex liquid phase at any temperature using the same model 
but with different geometry. 

In the disordered case, contrary to Kienappel and 
Moored, we showed that the glassy COCD phase ex- 
ists at low temperatures. In this phase the transla- 
tional invariance order parameter is found to decay as 
5(G,r) ~ e -»7G ln r as predicted theoretically by Car- 
pentier and Le Doussal 35 . However, due to small system 
sizes, the possibility that it decays as a simple power law 
is not excluded. The clarification of this point requires 
simulations for considerably larger system sizes. Our re- 
sults indicate, however, that quasi-long range order of 
correlations measured by exist at low temperatures. 

In the glassy phase the phase coherent vortex glass or- 
der parameter decays as Cvg{ y ) ~ exp(— r/i? c ), where 
the correlation length R c depends not only on the tem- 
perature but also on the disorder. The stronger the dis- 
order the larger is R c suggesting that the phase-coherent 
vortex glass ordering becomes more and more favourable 
as the disorder is enhanced. 

The outline of the paper is as follows. In Sec. II we 
present the LLL approach in the Landau gauge. The 
transition from the vortex liquid to the ordered phase of 
the clean system is studied in Sec. III. The nature of 
ordering in the same system but this time with disorder 
is discussed in Sec. IV. The exponential decay of the 
phase coherent vortex glass is also presented in this sec- 



tion. Finally, in the last section we summarize our main 
results. 



II. MODEL 

Our simulation is based on the phenomenological 
Ginzburg-Landau model in the approximation of a uni- 
form magnetic field B. This is a reasonable approxima- 
tion, since the effective screening length A e // = 2A|/s 
diverges for vanishing film thickness s. Here Xl denotes 
the bulk screening length. The quenched disorder is in- 
troduced via the fluctuations of the local phase transition 
temperature. The Ginzburg-Landau free energy of the 
superconductor is given by 



F 



F c i + Fdi 



F rl = I d 2 r [a(T)\V\ 2 



2 1 1 



— |D*| 2 ] 

2m 



F d 



d 2 r a <5T c (r)|*(r)| ; 



(9) 



Here VP is the complex order parameter representing 
the macroscopic wave function of the superconducting 
electrons. D denotes the gauge invariant derivative 
D= —ih\7 — with A being the vector potential. 
B = V x A, e and m are the charge and mass of the 
electron, respectively. In the simulation we will use the 
the Landau gauge A = B(0,x, 0). We may go back to 
the symmetric gauge A = -§(— y, x, 0) by a simple gauge 
transformation W L = ^ s e lxy/l2 . a(T) = a (T - T c ) and 
(3 is a constant; a , /? > 0. T c denotes the zero-field mean- 
field transition temperature. 8T c (r) is real and Gaussian 
distributed with 



< ST c (r) > 
< ST c (r)ST c (r') > 



0, 



(10) 



<5{ (x) is a S~ function of width of the order of the zero 
temperature correlation length £o = ^/ \/2maoT c . The 
typical fluctuations ST c (r) of the mean-field transition 
temperature are then of the order ST C = CJT C . It is con- 
venient to collocate here the main characteristics of the 
model 
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Here -B C 2(0) denotes the upper critical field at T = 
and I the magnetic length. The density and the lattice 
constant of the flux line lattice are given by \ and oa, 

"u 

respectively. Gi is the Ginzburg number at T = T c . 

In the LLL approximation one expands the order pa- 
rameter ^ in terms of eigenfunctions of the operator 
a + D 2 /2m of the lowest Landau level only. This restric- 
tion is a reasonable approximation provided that fluctu- 
ations in higher Landau level channels can be neglected. 
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The precise range of applicability of the LLL approxima- 
tion is still under debat o 20 i 42 . However, one can expect 
that this approximation is valid for B > £? c2 /13^. Since 
the basic length scale of the LLL approximation is given 
by the magnetic length I = ^Jhc/2eB it is convenient to 
measure from now on all lengths in units of I, i.e. we 
replace everywhere f by x and j by y. 

We apply quasi-periodic boundary conditions to the 
order parameter inside a finite system with lengths L x 
and L y parallel to the x— and y— direction, respec- 
tively. V(y + L y ,x) = $(y,x) and V(y,x + L x ) = 
e'Kp(—iL x y/l)' l Sf(y, x). Such boundary conditions are nec- 
essary to result in a phase change by 2i:N l p when orbiting 



the system. N$ 



2tt1 2 



denotes the number of vortices 



in the system. We will denote L x /l = 3 1 ' 4> n 1 ' 2 N x and 
Ly/l = 27r 1 / 2 iV a /3 1/4 , where N x and N y are the num- 
bers of vortex columns and rows, respectively. The order 
parameter can then be written in the form 1 ^ 



*(r) = 



\ob\1^ z 

~jLy 



E E c, 

j — l S — ~ OQ 



(12) 



Using the expansion 1)12(1 we can rewrite F c i as follows 



In the last equation we have introduced the following 
notations 



olb = ot(T) + HeB /to = cxqT c (1 — b — t), 

t = —,X jiS = j2nl/L y + sL x /l (13) 



In the mean field theory the superconducting transi- 
tion of clean systems occurs at T C (B) — T c — fieBja^mc 
which is defined from the condition as — 0. Xj lS is the 
center of the cyclotron motion. The number of vortices 
Ncj, must be chosen to be an integer, the coefficients Cj 
are in general complex. 



= sign(a B )J2\ C i 



3 i/4 



OO OO 

E E 



+f(n s + 1/2) 



exp[ 



exp[ 



V3n(n d - 1/2) 2 



n B +n d ]C-[n s -n d 



]C[n s +n d ]C[ ns - nd + i] 



(14) 



where Cr„ 1 = C, 



jnod^risjNtp) ■ 



The only model dependent parameter of the model in the clean limit is e given by 



olb 



1 (b + t-1) 

01/3TV2 ~ 2Gii bVH 1 / 2 ' 



(15) 



which measures the distance to the mean-field phase boundary e = 0. Here b denotes the magnetic field in units of 
the mean-field upper critical field B C 2 at T = (compare l|ll[l). Below we will find the phase transitions at a critical 
value e c of e which can be inverted into a band of transition lines b c (t) parameterized by the Ginzburg number of the 
corresponding superconductor: 



b c (t;Gi,(= 0) = b c {t,e c Gii) = l-t + x 2 c t- x c y/2t(l - t) + >* 2 i 2 , >c c = V2e c Gii . 



(16) 



Finally, 



f(n s ) = Erf(N x 3 1 / 4 (27T) 1 / 2 - n s ^ 1 3 1 / 4 (2 7 r) 1 / 2 ) + Sr/(n s iV- 1 3 1 / 4 (27r) 1 / 2 ), 



(17) 



where Erf(x) is an error function. It should be noted that different groups used different notations for the dimen- 
sionless parameter e. Our e = t (Ref. GJ) = t/2 (Ref. Gjj) = a T /V2 (Ref. E3) = g/V^ (Ref. E3). 

To express the disorder term Fdi S in terms of the LLL coefficients Cj , we expand the random Gaussian disorder in 
renormalized Hermite polynomials Uk(x) and harmonics 



<5T c (r) 



OO OO 



(Lyl) 



171 E E a km e^ imyl ' L «u k {x), u k {x) = (2 fe fcV^T 1/2 e- K ' 2 H k {x). 



(18) 



k=0 m=—oo 



Here H k (x) are Hermite polynomials and a km are complex random numbers. To satisfy ST c (r)* = ST c (r) one has to 
choose akm such that a kl - m = a kim . Using the standard orthogonality relations one can show that 5T c (r) satisfies 
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Eq. (fTU|) if < a km >= and < a k 



>= 



,.<6 



kk' 



From Eqs. l(T2*)l and (fTHf) we have 



fc=0 ra s ,n d 



ak,2nd-lC[ na+nd ]C* ns _ nd + 1 }L(k, Tl s + 
I 



1 



(19) 



Here 

L(k,n s ,n d ) 



^3 1 /2„2 



Jo 



dxuk{x)e 



-(a 



(20) 

and the dimensionless parameter £ which controls the 
relative disorder strength has the form 



c 



27r l/2 G jl/4^1/2 



7r 1 /2(i_ i _6)' 



(21) 



Substituting (|12[) into Eq. 0} one can express the 
Fourier transform of the density-density correlation func- 
tion in the following form (again we measure k~ l in units 
of I, i.e. we replace kl by k) 



where 



-k 2 /2 



< |A(k)P > , (22) 



A ( k ) =Y. C h-n v+ N^e 



L (n y -2j) 



(23) 



In what follow we will use < |A(k)| 2 > to characterize 
the Bragg peaks. 

Our MC simulations are performed for coefficients 
Cj G C** which are updated by the standard Metropo- 
lis algorithm. The candicate for the new configuration 
is generated by Cj — > Cj + eAC where AC is a com- 
plex number which is randomly chosen from the region 
|ReAC| < 1 and |ImAC| < 1 in the complex plane, e 
is chosen to be of order 0.1 so that the acceptance ratio 
is 0.3 ~ 0.5. The physical quantities are measured every 
20 - 50 MC steps. 



III. CLEAN SYSTEM 

The aim of this section is twofold. First, we want to 
check our code for the clean system (£ = 0) which was 
studied in detailiS*^. Second, we consider the spatial be- 
havior of the translational invariance correlation function 
which have not been studied previously by simulations. 

The simulations were carried out for systems of five dif- 
ferent sizes, each containing 4 2 ,6 2 ,8 2 ,10 2 and 12 2 num- 
ber of vortices. It took 3 x 10 4 ~ 10 5 MC steps to reach 
the thermal equilibration. The physical quantities were 
calculated over 10 5 ~5x 10 5 MC steps. 




FIG. 1: The wave vector dependence of the structure factor 
for the clean system with f = 0. iV^ = 100, e 2 = 20 and 100. 
Ox is measured in units -k 



Fig. n shows the temperature dependence of < 
|A(&)| 2 >=< lA^O)! 2 > for = 100, e 2 = 20 and 
100. At low temperatures (large — e) this quantity has 
the sharp peaks which are indicative for the existence of 
the quasi vortex lattice. At high temperatures the Bragg 
peaks are smeared out suggesting that we are in the vor- 
tex liquid phase. 

In order to obtain the quasi solid-liquid transition 
temperature we monitor the scaling of the renormalized 
structure factor at the maximum position q x = G = 

(4tt)" '' 



3l/4 , < | A(G?)| >, which is defined as follows 

< |A(G)| 2 > = Nj, < |A(G)| 2 > / < |A(0)| 2 > . (24) 

Fig. [SJshows the temperature dependence of < |A(G)| 2 > 
for various values of N$. The transition temperature is 
defined from the point where the curves splay out and 
we find e 2 = 43 ± 2. This can be compared with previous 
estimates e 2 = 43.5 ± 1.0 by Hu and MacDonald 1 ^, e 2 ~ 
50 by Tesanovic and XingH, and e 2 « 49 by Kato and 
Nagaosa 1 ^. The discrepancy between estimates of various 
groups is probably related to the fact that the scaling 
regime was not reached due to small system sizes used 
up to now. Since this regime, according to Kato and 
Nagaosa 19 , corresponds to > 36 2 and the CPU time 
grows as jVJ it is beyond our computational facilities to 
reach it. 
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FIG. 2: The temperature and size dependence of the structure 
factor for the clean system. = 16, 36, 64, 100 and 144. 
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FIG. 3: The temperature-field phase diagram for the clean 
system. The dashed line denotes the mean field boundary 
between the vortex liquid and the quasi-solid phases. The 
dot-dashed, solid and dotted lines correspond to the phase 
boundary for Ginzburg number Gi = 10 -7 , 10 -2 and 1, re- 
spectively. 



FIG. 4: The distance dependence of the translational invari- 
ance order parameter S(G, r) for the clean system. Closed 
circles, hexagons, squares and triangles correspond to e 2 = 
10, 50, 70 and 100, respectively, r is measured in units of I. 



o 

V 



O 
<f 

V 

Z 




FIG. 5: Dependence of the structure factor on N<f, for the 
clean system. Open triangles, squares, and closed triangles, 
squares and hexagons correspond to e 2 = 10, 30, 43, 70 and 
100, respectively. Dotted lines are linear fits and the solid 
line corresponds to rjc = ~ from the KTBHNY theory. 



Fig. shows the phase boundary between the vortex 
liquid and the quasi-solid phase using Eq. (|16|) . We 
have chosen the Ginzburg number Gi = 10~ 7 which is 
typical for low-T c superconductors, and Gi = 1CP 2 and 
1 for high-T c materials^. Since the thermal fluctuations 
enhance with Gi, the quasi-long range ordered vortex 
lattice shrinks as one increases the Ginzburg number. 

The distance dependence of the translational invari- 
ance order parameter S(G,r) for = 144 and for var- 
ious values of e is shown in Fig. At low temperatures 
we find a power law behavior of the correlation function 
for translational long range order 

S(G,t) ~ r-" G ( T ) (25) 



whereas at high temperatures an exponential dependence 
on r occurs. 

One can define t]g from the distance dependence of 
5(G,r) shown in Fig. 0] However, in order to mini- 
mize finite size effects we calculate from the depen- 
dence of < |A(G)| 2 > on N,p in the quasi solid phase: 

< |A(G)| 2 >~ n^ vg/2 (see Ref. M). In the liquid 
phase < |A(G)| 2 > docs not depend on A^. Fig. |S] 
plots In < |A(G)| 2 > versus hxZV^. The dependence of 
t]g(T) on parameter e is presented in Fig. At the 
melting temperature rjci(T m ) = 0.31 ± 0.03 which co- 
incides with the value of ^ from the KTBHNY theory. 
The fact that our results agree with the KTBHNY pre- 
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FIG. 6: Dependence of tjg(T) on e~ 2 for the pure system. 
The arrow indicates the melting temperature which corre- 
sponds to e 2 . =43. 



dictions suggests that the first order transition is weak 
and the second phase transition theory still applies. The 
weakness of the first order transition was supported by 
analytical calculations of Hikami et who obtained a 
small difference between free energies of the liquid and 
quasi-long range lattice phases at the melting tempera- 
ture. Since in the LLL approximation screening effects 
are neglected, A — + oo and hence, as follows from Eq. 
(J5J, t)g = ~ jz- The results shown in Fig. demon- 
strate that, in agreement with the theoretical prediction^, 



IV. DISORDERED SYSTEM 

In this section we will study the effect of weak disorder 
on the flux lattice melting transition and on the behavior 
of different correlation functions. Fig. [7| shows a typical 
snapshots of vortex cores obtained for the clean case (£ = 
0) and two values of the disorder strength: £ = 0.03 and 
0.2 after 500000 MC steps. The runs were done at e 2 = 60 
for the system of = 64 vortices and are long enough 
so that the equilibrium was reached. For the uniform 
(C = 0) and weakly disordered (C = 0.03) cases we find 
the slightly distorted Abrikosov lattice. The dislocations 
are clearly seen in the case of £ = 0.2 and the quasi-solid 
lattice phase cease to be exist. Thus, in the weak disorder 
strength limit dislocations can be neglected. We restrict 
ourselves to this limit and concentrate on two values of 
disorder £ = 0.01 and 0.03. Some modest simulations 
have been performed also for £ = 0.02 and 0.04. 

The MC simulations were carried out for systems of 
= 4 2 ,6 2 ,8 2 10 2 and 12 2 vortices. Depending on N^, 
10 5 — 2 x 10 6 MC steps were used and half of them were 
spent on equilibration of the system. The disorder aver- 
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FIG. 7: Snapshots of the vortex positions for ( = 0, 0.03 and 
0.2 obtained at e 2 = 60 after 500000 MC steps. Light spots 
correspond to vortex cores. We choose N$ — 64. 
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FIG. 8: The wave vector dependence of the structure factor 



for the disordered system with £ = 0.01 .N$ = 100, e 
and 70. 



30 



age is done typically over 10 — 40 samples. 

Fig. |S| shows the temperature dependence of the in- 
tensity of Bragg peaks defined by < | ^ ( fZtc ) 1 2 > for high 
(e 2 = 30) and low (e 2 = 70) temperatures, £ = 0.01 and 
N(j, — 100. As in the clean case, at low T's we still have 
sharp peaks characterizing the ordering of the vortices. 

The transition temperature to the vortex liquid phase 
may be defined from the finite size scaling analysis of 
the structure factor. Fig. [5] shows the temperature de- 
pendence of < |A(G)| 2 > given by Eq. I|24(l for vari- 
ous system sizes and ( — 0.01. The curves splay out 
at e 2 = 48 ± 2 indicating that the disorder reduces 
slightly the transition temperature to the liquid phase. 
For C = 0.03 we have e 2 = 57 ± 2 (the results are not 



shown). 
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FIG. 9: The temperature dependence of the structure factor 
for the disordered system with ( — 0.01. The curves splay 
out at e 2 = 48 ± 2. 
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To calculate from the critical value e c (Cj) ; 
the true phase boundary b c (t, Gi, Q for the four given 
values of Q (Ci = 0.01, £2 = 0.02,4 = 0.03 and 
(4 = 0.04) we use the relation (|21|) . Since the max- 
imal value of t is one, the maximum value of £ is 
given by Q(Gi) = 2y/7rGiiQe c j. Having chosen for a 
given Ginzburg number Gi the values of £ smaller than 
(j(Gi), we find that the relation (|21(l defines a tem- 



perature tj 



b c {t ,e cj Gi*). 



from which we get 



Here function b r is defined 



b c (tj,GiX) 
by Eq. ftty. 

Fig. HOIshows the b—t phase diagram for the disordered 
case for Gi = 0.01. The true disorder strength is equal 
to £ = 0.06. The dashed line corresponds to the clean 
system. Obviously, the fluctuations due to disorder shift 
the transition line to lower magnetic fields. The inset 
shows the dependence of the Ginzburg-Landau energy 
distribution function P{E) on E for N$ — 144 and e 2 = 
50. The double peak structure gives strong evidence that 
the transition is first order. 

The question we ask now is what is the critical disorder 
strength ( c above which the COCD glass phase disappear 
at any nonzero temperature. To answer this question one 
has to, in principle, find the dependence of the effective 
critical temperature e 2 on £. But for large values of £ it 
is very difficult to equilibrate the system at low temper- 
atures and one cannot locate e 2 . Therefore, we restrict 
ourselves to a few values of the disorder strength and plot 
1/e 2 versus £ as shown in Fig. ^] The interpolation to 
1/e 2 . = which corresponds to the real zero temperature, 
gives Cc ~ 0.11. 

Fig. HI shows the ]xlS(G,t) - lnr plot for S(G,r) 
defined in © of the N$ = 12x12 disordered system with 
C = 0.01. For e 2 > e 2 = 48 lnS*(G, r) decays linearly, 
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FIG. 10: The t — b phase diagram for the disordered system 
with disorder strength £ = 0.06. We chose Gi — 10~ 2 . The 
closed squares denote the position of the melting line in the 
presence of disorder and the dotted line represents the case 
without disorder. In the disordered case the glassy COCD 
phase exists at low temperatures. The inset shows the dis- 
tribution P(E) of the Ginzburg-Landau energy measured in 
units of the mean-field energy Emf = N^ksTe 2 / (3a, where 
Pa is the mean field value of the Abrikosov ratio. P(E) was 
obtained for = 144 and e 2 = 50 
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FIG. 11: The dependence of e 2 on C The interpolation by a 
straight line gives £ c ~ 0.11. 



except for e 2 = 50 where the decay is much faster. A 
slightly better fit can be obtained if one plots \nS(G, r) 
versus ln 2 r as shown in Fig. ^] In this case S^G,]*) ~ 
exp(— t)q, In 2 7'), in agreement with the prediction of the 
COCD theorjiS^. The exponent »7g(^) depends not 
only on T but, as we will see below, also on the disorder 
strength . 

Fig. 1141 shows the dependence of the translational in- 
variance order parameter S{G, r) on r for Q = 0.03. In 



this case we still have the roughness ~ln r but the decay 
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FIG. 12: The distance dependence of the translational invari- 
ance order parameter S(G, r) for the disordered system with 
iV = 144 vortices. £ = 0.01, closed triangles, squares, open 
hexagons and closed circles correspond to e 2 = 30, 50, 60 and 
70, respectively. 
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FIG. 13: The same as in Fig. 1121 but data are plotted versus 
ln 2 r. 



is faster than in the £ = 0.01 case. 

Fig. E| shows the temperature and disorder depen- 
dence of exponent fjc(T) for N$ — 144. From the simula- 
tion we find that fjc increases with the disorder strength. 
The COCD theory, on the other hand, predicts a univer- 
sal disorder independent value of ?7g, provided one is in 
the asymptotic region where one sees true fixed point 
behavior. We interpret our numerical result as a cross- 
over effect: because of the relatively small system sizes 
one is probably not yet in the asymptotic region of the 
COCD theory (L s < L co ) and hence an increase of the 
disorder will results in a faster decay of the correlation. 
The other observation is that fjc increases with increas- 
ing temperature. The same remains true for tjg if we fit 
S(G,r) ~ e-^inr ( compare Fig £2} The result for 



FIG. 14: The same as in Fig. but for < = 0.03. Closed 
triangles, open hexagons, closed circles and open squares cor- 
respond to e 2 = 30, 60, 70 and 80, respectively. 



?7g is opposite to the expectation from the COCD the- 
ory in the asymptotic regime. However, our result can 
be understood again if we assume that the main effect 
for the roughness of the vortex lattice still comes from 
thermal fluctuations. They alone lead indeed to the in- 
crease of fjc with temperature. Thus, we interpret the 
result of Fig. ^| as a combined effect of thermal fluctua- 
tions and disorder on scales L s < L co . The present data 
does not allow for a clear distinction between an lnr and 
ln 2 r-behavior of ln5(G, r). 




FIG. 15: The temperature and disorder dependence of ?Ja(T) 
characterizing the decay of the translational invariance corre- 
lation function in the super-rough phase. N$ = 144, ( = 0.01 
(closed squares) and f = 0.03 (open hexagons). Error bars 
are mainly due to sample-to-sample fluctuations. 

We next study the positional vortex glass ordering 
characterized by the order parameter Spg (G,r) © 
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which we try to fit as 

In £pg(G, r) = —t)pg In r + const. 



(26) 



where rjpc is an exponent. The COCD theory predicts an 
lnr-behavior— with t\pq ~ T. Fig. 1 161 shows the distance 
dependence of Spg(G, r) versus lnr for ( = 0.01 and 0.03. 
At low temperatures the fit given by Eq. (|26|l works 
well and the corresponding exponent rjpc is presented 
in Fig. El As the temperature is lowered the tendency 
to the ordering gets more and more enhanced and rjpc 
decreases. Fig. 1 171 shows that t]pq is not sensitive to the 
disorder strength. Both findings are in agreement with 
the theoretical prediction 44 suggesting a quasi-long range 
glass ordering of e lGu . 
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FIG. 17: Temperature dependence of t/pg for £ = 0.01 (closed 
squares) and 0.03 (open circles). 



In r 



FIG. 16: LaS P a(G,r) -lnr plot for = 144, f = 0.01 (up- 
per pane) and £ = 0.03 (lower panel). For £ = 0.01 closed tri- 
angles, squares, open hexagons and closed circles correspond 
to e 2 = 30, 50, 60 and 70, respectively. For ( = 0.03 closed tri- 
angles, squares, open hexagons and closed circles correspond 
to e 2 = 30, 60, 70 and 80, respectively. 

Next we consider the vortex glass correlation function 
for a phase coherent vortex glass, Eq. |(HJ). In order to 
compute the phase-sensitive correlation function CVc( r ) 
we go to the symmetric gauge to obtain the phase of the 
order parameter. The typical spatial distribution of the 
phase in the symmetric gauge is shown in Fig. 1181 It 
has essentially the same shape as presented in the work 
of Brandt. 

We fixed the order parameter at the center of the rect- 
angular as shown in Fig. El as *(0) and compute CVc( r ) 
with the help of Eq. |JSJ, where |r| is the distance to this 
center. To improve statistics we average over 360 differ- 
ent directions of r keeping |r| = r fixed. 

Fig. 1 191 and I2U1 show the spatial behavior of the vortex 
glass order parameter Cvc{ r ) given by Eq. (JSJ for £ = 
0.01 and 0.03. In the low temperature region straight 
lines correspond to the fit 

C VG (r) ~ exp(-r/i? c ) , (27) 




FIG. 18: Space distribution of the phase of the order param- 
eter in the symmetric gauge. The darker region the larger 
is the value of the phase. The snapshot was obtained for 
= 16 after 3xl0 5 MC steps at e 2 = 50. The boundary 
between the dark and the bright region corresponds to the 
phase jump from ip = 2n to <p = 0. Going around the system 
the total phase change is 2irN<j > . 



where R c is the correlation length. At high temperatures 
Cvg{ v ) decays faster than exponential. 

Fig. 12 II shows the temperature dependence of R c . As 
the temperature decreases, the system gets more ordered 
and R c increases. The same tendency shows up if one 
increases the disorder strength. However we do not see 
any sign of a true or quasi-long range order of the phase 
coherent vortex glass correlation function. 
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FIG. 19: The spatial behavior of the vortex glass order pa- 
rameter CVg(i") in the \nCvair) — r plot. The closed trian- 
gles, squares, hexagons, open circles and squares correspond 
to e 2 = 10, 30, 50, 60 and 70, respectively. Straight lines are 
linear fits logCvc(R) ~ R for e 2 = 50,60 and 70. N = 144 
and C = 0.01. 
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FIG. 20: The same as in Fig. E2 but for £ = 0.03. The 
closed triangles, squares, hexagons, open circles and squares 
correspond to e 2 = 10, 30, 60, 70 and 80, respectively. The 
straight lines are fits for e 2 = 60, 70 and 80. 

V. CONCLUSION 

We have studied the phase diagram and different cor- 
relation functions in a 2D type II superconductor by 
Monte Carlo simulations using the LLL approximation. 
For the clean case, in accord with the previous results 
of other groupsi2iiii2iiS, we have shown that, the quasi- 



long range ordered vortex lattice phase exists at low tem- 
peratures. The exponent ijg was found to be propor- 
tional to e~ 2 as expected from the KTBHNY theory of 
the dislocation driven melting transition. This result sug- 
gests that the transition from the quasi-solid phase to the 
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FIG. 21: Temperature dependence of R c for the phase coher- 
ent vortex glass order parameter in the COCD glass phase. 
JV = 144, C = 0.01 and < = 0.03. 



liquid one is weakly first order. It is not excluded that 
the agreement between our results and the predictions of 
the KTBHNY theory is just an artifact of the finite size 
effect. 

In disordered systems there is a first order transition 
from vortex liquid phase to a COCD phase character- 
ized by the ln 2 r-behavior of the roughness. This phase is 
however only stable for small enough disorder strength, 
£ < C c ph 0.11. In the COCD phase the positional vor- 
tex glass correlation function decays as power law. In 
agreement with the COCD theory, the corresponding ex- 
ponent 7]pg grows with T linearly and remains almost 
unaffected by the disorder. 

The phase coherent vortex glass correlation function 
decays with r exponentially (see Eq. I|27[)) where the cor- 
relation length R c increases with disorder indicating that 
the disorder favours the vortex glass ordering. However 
the phase coherent vortex glass is suppressed by thermal 
fluctuations in weakly disordered 2D systems. 
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